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ABSTRACT 

The one-loop supersymmetric QCD quantum effects on the width of the un¬ 
conventional top quark decay mode t —»• H + b are evaluated within the MSSM. 
The study of this process is useful to hint at the supersymmetric nature of the 
charged Higgs emerging from that decay. Remarkably enough, recent calcula¬ 
tions of supersymmetric corrections to Z-boson observables have shown that 
the particular conditions by which the decay t —> H + b becomes competitive 
with the standard decay t —> W + b have a chance to be realized in nature. 
This further motivates us to focus our attention on the dynamics of t —> H + b 
as an excellent laboratory to unravel Supersymmetry at the quantum level in 
future experiments at Tevatron and at LHC. 





The recent discovery of the top quark at Tevatron B0 is reckoned to be the lat¬ 
est major event in the world of elementary particle physics. The weighted average of 
the CDF and .DO measurements yields a mass, mt = 180 ± 12 GeV, which is in good 
agreement with the Standard Model (SM) global analyses of the electroweak precision 
data. It can also be used to sharpen the prediction of the SM Higgs mass, Mh , in a 
range centered at ~ 100 GeV or below 0. Far from being a final confirmation of the 
SM, the finding of the top quark raises many questions on the nature of the spontaneous 
symmetry-breaking mechanism (SSB) that go beyond the SM, in particular whether the 
SSB is caused by fundamental scalars or by some dynamical mechanism postulating a 
new species of strongly-interacting fermions (such as in technicolour like models ||]), or 
perhaps involving the condensation of the top quark itself (such as in topcolour models 
II). In this paper we shall adhere to the extensions of the SM associated with elemen¬ 
tary scalars, specifically to the Minimal Supersymmetric Standard Model MSSM [|fl. The 
global fit analyses of precision data within the MSSM chart an increasing trend of compat¬ 
ibility with the MSSM [[?]]. They lead to a top quark mass, m t = 165 ± 10 GeV, consistent 
with the experimental measurements. Moreover, the richer structure of the top and Higgs 
sector of the MSSM gives rise to a very stimulating top/stop-Higgs/higgsino dynamics. 
The archetype example of it could be the non-standard decay of the top quark into a 
charged Higgs: t —► H + b. Although there are in principle many new exotic decays of 
the top quark in the MSSM, perhaps the latter mode is (if open) the closest one to the 
canonical mode t —> W + b in the SM and therefore the less difficult to handle from the 
experimental point of view. 

Apart from its obvious interest on its own, a further motivation || to consider the 
analysis of the charged Higgs decay of the top quark stems from the recent results of 
Z physics, particularly of the observed anomaly in the ratios Rb, R c as well as in the 
value of a s (Mz ) ||. These anomalies are discrepancies (at the 2 — 3a level) between 
the measured values of these observables and the corresponding predictions of the SM. As 


shown by early calculations of supersymmetric (SUSY) radiative corrections to Rb the 
theoretical prediction of this observable can be made in better agreement with experiment. 
More recently, there has been a flood of renewed interest in this subject (see e.g. m-\ ni) 
and it has been possible to establish in a more precise way the particular conditions by 
which the MSSM is able to cure, or at least to alleviate, the u Rb crisis”. The corresponding 
MSSM analysis of the ratio R c (in correlation with the ratio Rb) was first presented in 
Ref. [141. The upshot is that both the Rb and R c “crises” can be solved within the MSSM 
provided that tan/3 is large enough [] and there exists a supersymmetric pseudoscalar 

1 A possible solution also exists for tan (3 < 1, but it is not favoured by model building |l(j] and it is 
not relevant to the present analysis. 
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Higgs as well as some superpartner all of them in the 50 GeV range -in agreement with 
the most recent global fit analyses [0. In these conditions it turns out that one can 
simultaneously provide a supersymmetric explanation |T^J of the longstanding mismatch 
between the low-energy and high-energy determinations of a s (M z ) ||18|. All in all, these 
results bring in an independent incentive from the high precision world of Z physics to 
choose our SUSY parameters in the region of large tan/? and moderate charged Higgs 


mass. Indeed, from the well-known Higgs mass relations in the MSSM [^9] and assuming 
that a light CP-odd (“pseudoscalar”) Higgs mass exists in the 50 GeV ballpark^, 
it follows that there must be a charged Higgs companion of M#± — 100 GeV. In these 
circumstances the decay t —> H + b becomes competitive with the SM decay t —> W + b 
and it should be possible to identify it by tagging violations of lepton universality caused 
by the presence of an excess of final state r-leptons associated to the subsequent Higgs 
decay -as commented at the end of Ref. [IT - 


In view of the potential interest of the decay mode t —*■ H + b, one would naturally like 
to address the computation of the strong virtual corrections to its partial width. Of these, 
the conventional QCD corrections have already been considered in detail in Ref. [pOfl and 
they turn out to be sizeable and negative (of order —10%). Although they are blind to the 
nature of the underlying Higgs model, they need to be subtracted from the experimentally 
measured number in order to be be able to probe the existence of new sources of quantum 
effects beyond the SM. These effects may ultimately reveal whether the charged Higgs 
emerging from that decay is supersymmetric or not. 

In this paper we compute the strong SUSY radiative corrections to the partial width by 
paying special attention to the aforementioned privileged region of the MSSM parameter 
space. The analysis of the larger and more complex body of SUSY electroweak corrections, 
namely the corrections mediated by squarks, sleptons, chargino-neutralinos and the Higgs 
bosons themselves, will be presented elsewhere PQ To compute the one-loop QCD 
corrections to T t = T(t — > H + b) in the MSSM, we shall adopt the on-shell renormalization 
scheme where the fine structure constant, a, and the masses of the gauge bosons, fermions 
and scalars are the renormalized parameters: (a, M w , M z , M H , rrif, M S usy , •••) |25|. The 
interaction Lagrangian describing the tbH ± -ve rtex in the MSSM reads as follows: 

gV tb 


£Hbt ~ 


H b [m t cot /? Pft + rrii, tan (3 P L ] t + h.c., 


( 1 ) 


\/2M w 

where Pl,r = 1/2(1 =F 75 ) are the chiral projector operators, tan (3 is the ratio between 

2 It is noteworthy that at high tan/3 the approximate phenomenological lower limit W/jO 7 50 GeV 
on the light CP-even Higgs mass of the MSSM translates into m^o > 50 GeV for the CP-odd mass. 
Recent global fit analyses also favour a light Higgs mass in the SM and a light CP-even Higgs mass in 
the MSSM @. 

3 The study of the corresponding supersymmetric quantum corrections to the canonical decay t —> W + b 
has been presented by some of the authors in Refs.|2^| and J2|] (See also Ref. |2^1 ). 
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the vacuum expectation values of the two Higgs doublets of the MSSM 0 and V tb is the 
corresponding Kobayashi-Maskawa matrix element-hereafter we set V t b = 1 ( V t b = 0.999 
within ±0.1%, from unitarity of the KM-matrix and assuming 3 quark families). 

There are no oblique strong supersymmetric corrections at 1-loop order. The non¬ 
oblique vertex corrections originating from gluinos and squarks (stop and sbottom species) 
are depicted in Fig.l. The SUSY-QCD interaction Lagrangian relevant to our calculation 
is given, in four-component notation, by 


£ = 


9s 

V2 


Ql (Ar)ii 9 r Pl q ] ~ <f(K)ij Pl 9 r q R ] + h.c., 


( 2 ) 


where g r (r = 1,2,..., 8) are the Majorana gluino fields, (A r )ij(i,j = 1,2,3) are the Gell- 
Mann matrices, and q' a = {c/l, q R } are the weak-eigenstate squarks associated to the 
two chiral components PL, R q ; they are related to the corresponding mass-eigenstates 
% — {qii^} by a rotation 2x2 matrix (we neglect intergenerational mixing): 


q'a 


= 


id)- 

ab qb, 


R( q) = 


COS 8q 

- sin 6 n 


sin 0g 
cos 6 n 


(q — 


/ q owe wq j 

These rotation matrices diagonalize the corresponding stop and sbottom mass 


Mj = 


M\ — 


M~ ± ml + cos 2/3(1 - § s^) Mf 


m t M\ 


LR 


v m Mi R 

M~ ± ml ± cos 2/3 ( 1 


M~ r ± m 2 ± | cos 2/3 s'ly Mf. 



2 + 3 S W ) 


bL 


m b M b LR 


mb M\ 


b 

LR 


M- 2 ^ ± ml — | cos 2/3 s( v M. 


( 5 ) 


Z , 


with 


M i LR = A t — g cot (3 , M b LR = Ab — /i tan/3, (6) 

/i being the SUSY Higgs mass parameter in the superpotentialQ. The A tjb are the trilincar 
soft SUSY-breaking parameters and the Mq LR are soft SUSY-breaking masses 0. By 
SU (2)^-gauge invariance we must have Mf L = M^, whereas M^, M~ bR are in general 
independent parameters. Finally, we also need the interaction Lagrangian involving the 
charged Higgs and the stop and sbottom squarks 


P ml ~ ~ H (g LL b* L t L + g RR b* R t R + 9lr b R t L + 9rl b* L Hr) ± h.c., (7) 

where 


gLL = Mf, sin 2/3 — (ml cot (3 ± m 2 tan (3 ), 

9rr = —m t m b (ta,n (3 ± cot (3 ), 
g L R = -m b (g + A b tan/3), 

9rl = -m t (/i +A t cot (3). ( 8 ) 

4 Its sign is relevant in the numerical analysis. We fix it as in eq.(3) of Ref. p2|. 
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The one-loop renormalized vertex, A, is derived from the renormalized Lagrangian pins 
counterterms, C —> C + 8C, following the standard procedure |25[]. It can be parametrized 
in terms of two form factors Fr , Fr and the corresponding mass and wave-function renor¬ 
malization constants Srrif, 8Z[ R associated to the external quarks, viz. 

ig 


A = 


x/2M, 


w 


[m t cot f3 (1 + A R ) P r + m b tan/? (1 + A L ) P L \ 


(9) 


with 


A r = F R + - + -6Z b L + -6Z t R , 

H H mt 2 L 2 R ’ 

a l = u + — + Fz‘ L +Uz b R . 

m b 2 2 


In the on-shell scheme we haveQ 


8m n 


m n 


siK) + s?iK) , 2 , 

-- - + £«(">,) 


and 


6Zl„ = E»(m5)+mJ[E»'(mJ) + Ef i '(mJ) + 2E|'(m^)] 


( 10 ) 


( 11 ) 


( 12 ) 


In these equations we have decomposed the (real part of the) quark self-energy according 
to 

J2 f (p) = E{(p 2 ) i>P L + Z f R (p 2 ) i>P R + m f J2 f s (p 2 ) , (13) 

and used the notation E '(p) = dJ2(p)/dp 2 . 

From the renormalized amplitude ([|), the width T = F(t —> H + b) including the 
one-loop SUSY-QCD corrections is the following: 


r = r 0 1 + 


{l + Yp&fAJl + ^p&tA^ + YspflefA. + AH)]} , (14) 


N L v /J N 
where the corresponding lowest-order result is 


r 0 =(-^-A 1/2 (i,^,^%). 


\8ttV 2 m t 


2 5 2 
Wit m t 


We have defined 


A 1/2 (l, x 2 , y 2 ) = til- (i + !/) 2 ]( 1 - (X ~ s) 2 ] 


(15) 

(16) 


and 


N = ( m 2 + ml — M r+ ) ( m 2 cot 2 f3 + m 2 tan 2 (3) + 4m 2 m 2 , 
-/Vl = (m 2 + ml — M r+ ) ml tan 2 /3 , 

IV R = (ml + ml — M r+ ) ml cot 2 (3 , 

N lr = 2 m 2 t m 2 


b ■ 


(17) 


3 Our sign conventions for the self-energy functions are those of RefJ 
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Notice that, in contradistinction to electroweak one-loop calculations ||22||, an additional 
correction term A r does not appear in (|14|) due to the absence of one-loop SUSY-QCD 


corrections in fi- decay. 

The explicit contribution to the form factors F l ,F r from the vertex diagram of Fig.l 
is given by 

F l = 8Tra s iC F G ; h J R$*R { £(Cn - c 12 )m t + R$> R®C 12 m b + R%i*R%>C 0 mj ], 

777/^ tclll 

Fr = 8ira s iC F ° ab g (C u - C 12 )m t + R$*R%>C 12 m b + R®*R$>C 0 m 5 ]. 

m t co l fj 

(18) 

Here C F = (Nq — l)/2Nc = 4/3 is a colour factor. We have furthermore defined: 

G ab = RSR^gu, + R 2b R 2 I*9rr + R$RS*9lr + 9rl ■ (19) 


The three-point function notation is as in Ref. 22, 27], with the following arguments: 


C = C(p,p',v. 


( 20 ) 


As for the self-energies, 


El(p 2 ) = 8na s iC F \R { ?}\\B 0 -B l ), 

E^(p 2 ) = 8na s iC F \R^\ 2 (B 0 ^B l ), 

S|(p 2 ) = -8na s iC F ^-Re(R^ R^*)B 0} (21) 

m q 

where the two-point functions -defined also as in Ref. |^2|- have the following arguments: 


B = B(p 2 ,m| a ,m?). 


( 22 ) 


(A summation over squark indices is understood in eqs.([l^) and (^).) It is easy to 
convince oneself that the form factors (|T^) are to be HR-finite in SUSY-QCD, as indeed 
they are. The remaining contributions to the renormalized amplitude @ are immediately 
seen to cancel UV -divergences each other out. 

The numerical analysis of the strong supersymmetric corrections to T(t —> H + b) is 
exhibited in Figs.2-5. We present the results both in terms of the corrected Y and in 
terms of the relative correction with respect to the tree-level width, i.e. 


h = • ( 23 ) 

1 0 

The free parameters at our disposal lie in the mass matrices (f|)-(j§j)- Among the sfermion 
mass matrices, the stop mass matrix is the only one where a non-diagonal structure caused 
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by a sizeable mixing term is most likely to arise. For M\ jR = 0 this matrix is trivial, but 
since m t is large a nonvanishing M\ R naturally leads to a light mass eigenvalue, denoted 
by rnj A , whereas the other eigenvalue, rrif 2 , can be much heavier. As a matter of fact 
a light stop with a mass = O^Adz/Z) is still phenomenologically allowed [^6] 


In 


contrast, the off-diagonal element of the sbottom mass matrix ([]), being proportional to 
m b ~ 4.5 GeV, is expected to be small (unless M RR is very large). We shall treat the 
sbottom mass matrix in the simplest possible way compatible with the phenomenological 
bounds on squark masses [^6|- only scaped perhaps by the lightest stop. For definitiveness, 
unless stated otherwise, we shall assume that M RR = 0 (equivalently A b = /rtan/3) and 
that the two mass eigenvalues are equal (mg = m~ b2 = m b ) and constrained to satisfy 
m b > 150 GeV. As for A t , it will be treated either as an input parameter or it will be 
fixed once we are given Ad RR , /i and tan/A We remind that M RR is expected to preserve 
the inequality 

Mlr — 3 rn b L > ( 24 ) 

which roughly corresponds to a necessary, though not sufficient, condition to avoid colour¬ 
breaking vacua [^]. In the conditions described above, once m b , /i and tan/? are fixed, 
the stop mass matrix depends on only 2 parameters, e.g. (A t ,Mj R ), {Ad\ jR , rrqj, etc. For 
the strong coupling constant we used the value 


( 7 2 

= o.ll 

47T 


(25) 


which remains essentially constant within the CDF-DO ranges mentioned above. When¬ 
ever m t needs to be fixed, we take the central value m t = 180 GeV. 

A crucial parameter to be explored in our analysis is tan/?. In Fig.2a we plot the 
SUSY-QCD corrected Y = Y(t —> H + b ), eq.([T|), versus tan/? for n = +100 GeV and 
li = —100 GeV, and for given values of the other parameters. We see that T is very 
sensitive to tan/3 and that, barring the narrow interval tan/? < 1, the process t —> H + b 
becomes steadily competitive with the standard process, t —> W + b, in the large tan /? 
region, i.e. when tan/? is of the order of 30 — 40 ~ m t /'m b . As mentioned in the beginning, 
this is precisely the range singled out by the Z boson observables (T4|]. We also see from 
Fig.2a that, for tan/? > 30, Y starts to deviate from To quite manifestly. Therefore we 
shall choose tan /? = 30 as a representative value in the other plots. 

Highly remarkable is also the incidence of the parameter /r both of its value and of its 
sign. In fact, the sign of 5g happens to be opposite to the sign of /i and the respective 
corrections for fi and for —// take on approximately the same absolute value. In Fig.2b we 
deliver the correction itself, Sg, for /i = —100 GeH and for different values of the squark 
and gluino masses. The sign dependence of 5g suggests that two extreme scenarios could 
take place with the SUSY-QCD corrections to T(t —> H + b)\ namely, they could either 
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significantly enhance the, negative, conventional QCD corrections |20|, or on the contrary 
they conld counterbalance them and even result in opposite sign. 

In Figs.3a-3b we display Sg as a function of m,g and m*, respectively, for three values 
of the sbottom masses. We learn from Fig.3a that light gluinos of 0(1) GeV m yield, 
contrary to naive expectations, a rather small correction as compared to heavy gluinos of 
0 ( 100 ) GeV. It should be clear that these corrections do eventually decouple -as we have 
checked- for larger and larger gluino masses. Notwithstanding, the decoupling rate of the 
gluinos is particularly noticeable, for it happens to be so slow (Fig.3a) that it fakes for a 
while a non-decoupling behaviour. This trait is caused by the presence of a long sustained 
local maximum (or minimum, depending on the sign of /i) spreading over a wide range 
of heavy gluino masses centered at ~ 300 GeV. For this reason, heavy gluinos are in the 
present instance preferred to light gluinos. In Fig.4, we test the sensitivity of Sg to M H ±; 
it turns out to be very small, except near the uninteresting vicinity of the phase space 
border where r 0 is about to vanish. 

In all previous figures we have fixed A t = 0, 9 t = 7t/ 4 and m~ b < 150 Geld, so that 
the lightest stop mass was always > 190 Geld. In Fig.5 we relax these conditions and 
plot contour isolines Sg = const, in the rn^ )-plane, where A t and 6 t are variable. In 

particular, we approach the region of the lightest possible stop masses compatible with the 
strict LEP phenomenological bound. As expected, the corrections are larger the smaller 
is . Finally, let us mention that a non-vanishing mixing in the sbottom mass matrix 
does not alter at all the typical size of the corrections obtained here. It was only to avoid 
much cluttering of free parameters that we have treated so far that matrix in the most 
simple-minded form (M b LR = 0 and m bi = ). For example, if the mass-eigenvalues 

are mj; i2 = 200,250 Geld and fi = —100 Geld, tan/? = 30, then one has A b = —500 Geld 
and the corresponding correction reads Sg = +30%. In the same conditions, but choosing 
A b = 0, g is determined to be —83 Geld and Sg = +25%, etc. 

Some words on previous work are in order. We chart significant differences in our 
complete analysis as compared to preliminary calculations in the literature. In Ref.[|30| a 
first study of the SUSY-QCD corrections to t —> H + b was presented, but they neglect the 
bottom quark Yukawa coupling and as a consequence they are incorrectly sensitive to the 
high tan f3 effects. Indeed, the form factor which is associated to that coupling is by far 
the dominant piece of our numerical analysis in the large tan /? region. Furthermore, the 
impact from mixing effects and the incidence of the various parameter dependences were 
completely missed and only the simplest situation, characterized by degenerate masses, 
was considered^. The study of Ref. |R] also neglects the bottom quark Yukawa coupling. 
Thus the purported large effects claimed in the large tan /? region are not correctly justi- 


6 Notice that the assumption of stop masses equal to sbottom masses is incompatible with eqs. (§-(!)• 













fied. Moreover, in the framework of these two references, the lowest-order width is fully 
proportional to cot/3; thus, in such a context, finding quantum effects increasing with 
tan/? is rather useless since they result in corrections to an uninteresting, vanishingly 
small, tree-level width. Finally, the latter reference also neglects the wave-function and 
mass renormalization contributions which, numerically, are of the same order of magni¬ 
tude as the vertex contributions in the intermediate and small tan/? region, ffowever, 
where it has been possible to force an overlapping, we have found numerical coincidence. 

In summary, the SUSY-QCD contributions to the partial width of t —> H + b could 
be quite large (several 10%); and what is more, these corrections apply to a decay mode 
which, in the region of the MSSM parameter space prompted by the high precision Z 


boson observables fll4| , [TTjl , has an appreciable branching ratio as compared to the standard 
decay t —> W + b. Furthermore, we have found that the impact of the gluino corrections 
on t —> H + b could occur in two opposite ways, either by reinforcing the conventional 
QCD corrections or on the contrary by severely cancelling them out-perhaps even to the 
extend of reversing their sign!. Most remarkable, the potential size of these effects stems 
not only from the strong interaction character of the SUSY-QCD corrections, but also 
from the high sensitivity of t —> H + b to the (weak-interaction) SSB parameter tan /?. 
At the end of the day we must conclude that t —> H + b could reveal itself as the ideal 
environment where to study the nature of the SSB mechanism. It could even be the right 
place where to target our long and unsuccessful search for large, and slowly decoupling, 
quantum supersymmetric effects. In this respect it should not be understated the fact 
that the typical size of our corrections is maintained even for sparticle masses well above 
the LEP 200 discovery range. Theses features are in stark contrast to the standard decay 
of the top quark, t —> W + b, whose SUSY-QCD corrections are largely insensitive to 
tan/?|23|]. Fortunately, the next generation of experiments at Tevatron and the future 
high precision experiments at LHC may well acquire the ability to test the kind of effects 

1331,131 


considered here (Cf. Refs. 


Thus, in favorable circumstances, we should be 


able to disentangle the potential supersymmetric nature of the charged Higgs decay of 
the top quark out of a measurement of the top quark width/] at a modest precision of 
~ 5 - 10%. 


7 Or related observables, such as e.g. the -double and single- top quark production cross-sections, or 
the differential distributions of the lepton final states in given exclusive decay channels, etc. 
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Figure Captions 

Fig-1 SUSY-QCD Feynman diagrams, up to one-loop order, correcting the partial 
width T(t —> H + b). Each one-loop diagram is summed over the mass-eigenstates of 
the stop, sbottom squarks ( b a , t b -,a,b= 1 , 2) and gluinos g r ; r = 1 , 2,..., 8. 

Fig.2 (a) SUSY-QCD corrected r(£ —► H + b) as a function of tan 3, for two opposite 
values of /x, compared to the corresponding tree-level width, r 0 . The framed set of 
inputs is common to Figs. 2,3. The horizontal (dashed) line marks T 0 (t W + b) = 
1.71 GeV -the tree-level width of the standard process ( m t = 180 GeV); its SUSY- 
QCD corrections are generally small and essentially independent of tan [3 |p3| ; (b) 
The relative correction 5g as a function of tan f3 for p = —100 GeU and three values 
of mr = rrin 


9' 


Fig.3 Dependence of 8g upon (a) trig, including the light gluino region, for the same 
squark masses and /i as in Fig.2b; (b) 5g as a function of rn t (within the CDF-D 0 
limits) and remaining parameters as in (a). 

Fig.4 5g as a function of M H ±. Rest of inputs as in Fig.3. 

Fig.5 Contour plots of 5g in the mq)-plane for ji = —100 and tan (5 = 30. 

The mixing angle 9 t and A t are variable and the other fixed parameters are as in 
Fig.2a. The shaded area is excluded by the triple condition > 0, > M z / 2 


and eq. (|24D . 
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